Astronomy & Astrophysics manuscript no. Mueller6062 


February 5, 2008 


(DOI: will be inserted by hand later) 





The multi-component field topology of sunspot penumbrae 
A diagnostic tool for spectropolarimetric measurements 

D.A.N. Miiller 1 , R. Schlichenmaier 2 , and G. Fritz 2 ' 3 and C. Beck 2 



O 
O 



European Space Agency, Research and Scientific Support Department, 
c/o NASA Goddard Space Flight Center, Mail Code 612.5, Greenbelt, MD 20771, USA 
[ e-mail: dmueller@esa.nascom.nasa.gov 

2 Kiepenheuer-Institut fur Sonnenphysik, Schoneckstr. 6, 79104 Freiburg, Germany 
e-mail: schlicheOkis . uni-freiburg . de , cbeckOkis . uni-freiburg . de 

3 Arnold Sommerfeld Center for Theoretical Physics, Theresienstr. 37, 80333 Miinchen, Germany 
e-mail: derfritz@gmx.net 



C/3 



(N 
(N 



Draft version 



> 



Abstract. Sunspot penumbrae harbor highly structured magnetic fields and flows. The moving flux tube model offers an ex- 

■ planation for several observed phenomena, e.g. the Evershed effect and bright penumbral grains. A wealth of information can 
£f*l | be extracted from spectropolarimetric observations. In order to deduce the structure of the magnetic field in sunspot penum- 

■ brae, detailed forward modeling is necessary. On the one hand, it gives insight into the sensitivity of various spectral lines to 
^\ ' different physical scenarios. On the other hand, it is a very useful tool to guide inversion techniques. In this work, we present a 

generalized 3D geometrical model that embeds an arbitrarily shaped flux tube in a stratified magnetized atmosphere. The new 
semi-analytical geometric model serves as a frontend for a polarized radiative transfer code. The advantage of this model is 
that it preserves the discontinuities of the physical parameters across the flux tube boundaries. This is important for the detailed 
. shape of the emerging Stokes Profiles and the resulting net circular polarization (NCP). (a) The inclination of downflows in the 

O |- outer penumbra must be shallower than approximately 15°. (b) Observing the limb-side NCP of sunspots in the Fel 1564.8 nm 

line offers a promising way to identify a reduced magnetic field strength in flow channels, (c) The choice of the background 
atmosphere can significantly influence the shape of the Stokes profiles, but does not change the global characteristics of the 
resulting NCP curves for the tested atmospheric models. 
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1. Introduction velocity gradients are a necessary and sufficient condition to 

produce a NCP. 

s_j ■ Spectral line profiles in the penumbra are characterized by line _ . , , . ,.,.„, , 

. , ,. .ci - /r. i ii The advent or modern polanmeters l ike th e Advanced 

_C5 shifts and line asymmetries of their intensity (Stokes I) as well Po i ar1meter f ASP kknrnanich et al Il l99i at the NSO 

. . i i- « /0 i s-^ tt i t r\ * -ii oLOK.es r oianiiieLer irtor, oivuiiiaiiicii ei ai. iyy / ) ai me i^ou 

as in polarized light (Stokes Q, U, and V). A suitable measure , , _, _, , . — ~ — — XT " . , 

... j. o 1 t, jzi-i Vacuum Tower Telescope in Sunspot, New Mexico, and the 

of the asymmetry of Stokes- V profiles is the area asymmetry _, . „ T . , , . /Pr , T _, |„ . 7^ — , 

, . . r ,,• ■ •■ Tenenfe Infrared Polarimeter (TIP, Martinez Pilletetal. 1999) 

or net circular polarization, NCP, of a spectral line, which we , , , . . T . _ — /T ,^ T TO r"" - 

, _ and the Polarimetric Littrow Spectrometer (POLIS, Beck et al. 



define as 

NCPs j V(A)dA, (1) 

SA 



Eool at the Vacuum Tower Telescope (VTT) on Tenerife has 

brought new challenges to observers as well as theoreticians. 

It is now possible to routinely record maps of the full Stokes 

vector for a whole sunspot. This facilitates the study of spectral 

, , „ . r , , , signatures of magnetic flux tubes in the penumbra by analysis 

where the interval of integration, oA, encompasses the whole ° ^ NCP 

line profile. This quantity has the advantage that it is indepen- ° ma P s e 

dent of the spatial resolution of the observations since a convo- So far ' ^ s P ectral synthesis has been limited to atmo- 
lution of a line profile with a Gaussian preserves its area. The s P herlc models based on snapshots of themoving tube model 
first measurement s of the net circular polarization in sunspots simulations by | Schlichenmaier et al. | 1199J and to configura- 
were reported by lllling et al.l Jl975l) . lAuer & Heaslevl dl978i) tlons ln a two-dimensional plane ("slab models"). Such_slab 
showed that these observations could be explained by assum- models have been used earlier by|Solanki & Montayon| jl993>, 

ISanchez Almeida et all i 19961) and lMartmez Pilled OOOfJ) . and 
led to the picture of the "uncombed" penumbra: Approximately 



ing macroscopic velocity fields. Furthermore, they proved that 
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background field that is at rest. The sharp gradient or even dis- 
continuity of the velocity of the plasma between the inside of 
the tube and the static surroundings result in the observed NCR 
While two-dimensional models reproduced the observed NCP 
and its center-to-limb variation, they are limited to the few ge- 
ometric cases where a flux tube is exactly aligned with the line- 
of-sight (LOS). Thus, they do not give any information about 

the NCP at diff erent positions within the penumbra. 

We showed ISchlichenmaier et all20o1lMuller et al.l2002l 
hereafter M02) that a wealth of additional information is con- 
tained in NCP maps, which can be retrieved with the help of 
a three-dimensional model. Discontinuities in the azimuth of 
the magnetic field (as measured in the observer's coordinate 
system) lead to a characteristic pattern of the NCP within the 
penumbra. Specifically, it was shown that the apparent symme- 
try within the spot is broken due to the following fact: The dis- 
continuities in the magnetic field azimuth are different for two 
locations that are symmetrically located in the spot with respect 
to the "line-of-symmetry", i.e. the line connecting disk-center 
and sunspot center (Fig.Q. These azimuth discontinui t ies ha ve 
been confirmed observationally bv lBellot Rubio et"al 1 il2004h . 




A0<O -s 



Fig. 1. Pictorial representation of the difference in azimuth, 
A0(t/O, between the magnetic field vector inside a flux tube 
(solid arrows) and that of the background field (dashed arrows). 
The tubes are nearly horizontal, while the field vector of the 
background has a steeper inclination, but is also in the plane 
spanned by the spot axis and the flux tube axis. 



The impact of this geometric effect on the NCP strongly 
varies between spectral lines at different wavelengths. A com- 
parison of the NCP of different spectral lines can thus be used 
to diagnose the magnetic field inclination and orientation of 
flow channels in the penumbra. As an example, we consider 
the two iron lines Fel 1564.8 nm and Fel 630.25 nm. Both 
lines are Zeeman triplets and have similar Lande factors of 
g=3 and g=2.5, respectively. It turns out that the infrared line 
(Fe 1 1564.8 nm) is far more susceptible to the symmetry break- 
ing than the visible line (for details, cf. M02). 



2. The VTUBE model 

In order to gain a better understanding of the different fac- 
tors that determine the NCP and its spatial variation within 
the penumbra, we have constructed a 3D geometric model 
(VTUBE) of a magnetic flux tube embedded in a back- 
ground atmosphere. This model s erves as the frontend for a 
radiative transfer code (PI AM AG. iGrossmann-Doerthll 1 994h . 
Combining the two, we can generate synthetic NCP maps for 
any desired axisymmetric magnetic field configuration and ar- 
bitrary flux tube properties. The model has been built to offer 
a high degree of versatility, e.g. the option to calculate sev- 
eral parallel rays along the line-of-sight that intersect the tube 
at different locations. One can then average over these rays to 
model observations of flux tubes at different spatial resolutions 
and for different magnetic filling factors of the atmosphere. 
Furthermore, we can also take into account radial variations of 
the physical properties of the flux tube. Doing so, one can e.g. 
model the interface between the flux tube and its surroundings. 

In the simplest case, the VTUBE model reduces to a gen- 
eral model for a two-component atmosphere, i.e. an atmosphere 
with a static background component and a flow component. 
Thus, it can also be used for flux tube models other than the 
moving tube model considered here. 



surface normal 




Fig. 2. Magnetic flux tubes in the penumbra. For a given he- 
liocentric angle 9, the azimuth <p and the inclination y of the 
magnetic field vary with the location of the flux tube in the 
sunspot, characterized by the spot angle 4>- 

Figure|5]shows the relevant angles that appear in the model. 
For a given heliocentric angle 9, the azimuth <p and the inclina- 
tion y of the magnetic field vary with the location of the flux 
tube in the sunspot, characterized by the spot angle ifr. For the 
calculations presented in this paper, no averaging over parallel 
rays was performed, and the physical quantities were kept con- 
stant along the flux tube radius (thin flux tube approximation). 

2.1. The background atmosphere 

Apart from the details of the magnetic flux tube, the VTUBE 
model allows the user to specify the background atmosphere. 
We assume that the thermodynamic part of the atmospheric 
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model only depends on the geometrical height, z- Thus, tem- 
perature, T, pressure, p, and density, p, are prescribed in tab- 
ulated form, with data obtained either from semi-empirical or 
theoretical models. The magnetic field structure is described 
independently from the thermal structure. This is done because 
we are interested in the effect of varying physical parameters 
of the flux tube, and there is so far no self-consistent 3D MHD 
model of flux tubes in the penumbra available. For atmospheric 
models for which the magnetic vector field has not been de- 
termined, one has to take care that the parametrization of the 
magnetic field is consistent with the thermodynamic structure. 

We use the same parametrization as in M02 to model an 
axisymmetric magnetic field whose strength decreases radially 
as B(r) = Bo/(l + (r/r») 2 ) and whose inclination decreases lin- 
early as y'(r) = 7rr/(2ry) . In these expressions, r denotes the 
radial distance from the spot center, r* is the radial distance at 
which the magnetic field has decreased by 50%, and y' indi- 
cates the radial distance at which the field becomes horizontal. 
This re presents an extension to the work o flBeckers & Schroted 
( 1969) wh ich tak es into account more recent observations (e.g. 
iTitle et alii 19931) . according to which the magnetic field of a 
sunspot is not horizontal at the spot boundary (r — r», where r» 
is the radius at which the field strength has dropped to half its 
maximal value), but shows a small inclination with respect to 
the surface. Adopting values of Bo = 2700 G, r, = 16 800 km, 
and ry = 17400km yields an inclination of 3.3° at r = r„ 
and a vertical gradient of the magnetic field in the spot cen- 
ter of dBJdz ~ -0.25 G/km, which agrees fairly well with the 
value of dB 7 /dz ~ —0.3 G/km as determined from the model of 
Jjahn & SchmidJl994 hereafter JS94) for r = 0. 

In this work, we experiment with three different back- 
ground atmospheres. The first one is the penumbral model of 
JS94 which wa s used as a background atmos phere for the MHD 
simulations of ISchlichenmaier et alJ Jl998l) and for the NCP 
ca lculations of M02. The se cond one is the penumbral model 
of lBellot Rubi o et al.l which is based on high-resolution 

intensity profiles of the non-magnetic Fel 557.6 nm line, ob- 
tained with the TESOS spectrometer at the Vacuum Tower 
Te lescope on Tenerife. The t hird one is th e quiet sun model 
of iHolweger & Mulled ( 1 1974 modified bv ISchleicheil ( 1 19761) 
to account for the chromospheric temperature rise. This model 
was chosen to study the effect of a slightly higher temperature 
in the line-forming region of the atmosphere. 

Figure |3 shows the physical parameters along the line-of- 
sight for the three different atmospheric models with a hori- 
zontal magnetic flux tube. The upper right panel shows the 
temperature of these three atmospheric models as a function 
of optical depth. Between log(r) = -2 and log(r) = 0.5, 
the quiet sun model is u p to 570 K hotter than the model of 
IBellot Rubio et alJ d20od) and up to 860 K hotter than the at- 
mosphere of JS94. 

Since spectra yield information about the atmospheric pa- 
rameters as a function of optical rather than geometrical depth 
and since the NCP is very sensitive to the exact location of 
the flow channel along the line-of-sight, the location of the 
flux tube has been adjusted to allow a comparison between 
the different models. The top of the tube was placed at a rel- 
ative height of z, = -50 km for the JS94 model, z, = +82 km 



for the penumbral model of IBellot Rubioe^al 1 l|2006l) and 
Zt = +78 km for the quiet sun model of TscMeichei dl976l) . so 
that the top of the flux tube is located at an optical depth of 
log(r) = -1 for a vertical line-of-sight (0 = 0°). The example 
shows a magnetic flux tube with a flux of T = 5 ■ 10 16 Mx in 
the outer penumbra (r — 12 000 km), observed at a heliocentric 
angle of 9 — 30° and a spot angle of if/ = 0° (i.e. center-side 
penumbra). The flow speed inside the tube is vo = lOkm/s 
(this speed is used for all calculations in this paper unless 
otherwise noted) which results in a line-of-sight velocity of 
vws = Vo ' sin(0) = 5km/s. The minus sign indicates that 
the flow is directed towards the observer. Since the flux tube 
is aligned with the line-of-sight, the azimuth angle is zero. 

2.2. Building a realization of the model 

Once the background atmosphere has been defined, the follow- 
ing steps are performed to assemble a realization of the penum- 
bral flux tube model: 

1 . The flux tube axis is defined by a polygon or spline func- 
tion. 

2. The physical properties of the flux tube are defined and then 
interpolated along the flux tube axis. 

3. The viewing angle is defined by specifying the heliocentric 
angle, 6, and the location of the line-of-sight. 

4. The spectral line profiles are calculated. 

These steps are described in detail in the appendix. 

3. From Stokes-V profiles to NCP curves 

To illustrate how the shape of Stokes-V profiles is connected 
to a NCP curve, Fig.0]shows Stokes-V profiles and the corre- 
sponding NCP for a straight horizontal magnetic flux tube. The 
tube contains a magnetic flux of f = 5 ■ 10 1 6 Mx, is embed ded 
into th e penumbral atmospheric model of IBellot Rubi o et alJ 
(2006) and is observed at a heliocentric angle of 6 = 30°. The 
stack plot on the left side shows Stokes-V profiles for differ- 
ent positions along a radial cut at r = 12 000 km through the 
outer penumbra, specified by the spot angle, if/. The spot an- 
gle runs counter-clockwise from the line connecting the center 
of the spot and the center of the solar disk (see Fig. |2). For 
example, if/ — 0° corresponds to a position in the center-side 
penumbra where the Evershed flow is directed towards the ob- 
server. This can be seen by the blue-shifted component of the 
Stokes-V profile. The right plot shows the corresponding NCP, 
i.e. the integral over the Stokes-V profiles. As expected, it is 
seen that the NCP vanishes for if/ = 90° and if/ = 270°, where 
the line-of-sight is perpendicular to the flow. 

The Zeeman splitting of a spectral line is proportional to 
the square of its wavelength, A, while the Doppler shift in- 
creases only linearly with A. Due to this elementary difference, 
the effect of a flow channel in the line-of-sight of a magne- 
tized atmosphere is most clearly seen for longer wavelengths, 
e.g. in the infrared. In the upper and lower profiles in Fig. |4] 
(if/ € [330°, 30°]), the flux tube component of the Stokes-V pro- 
file is clearly visible as a small hump (or "line satellite") in the 
wing of the background component. In this case, the observer 
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Fig. 3. Physical parameters along the line-of-sight for three different atmospheric models. The geometrical location of the flux 
tubes has been adjusted so that the top of the flux tube is located at log(r) = -1 for a vertical line-of-sight (8 = 0°) for all three 
models. The example shows a magnetic flux tube with a flux of f — 5 ■ 10 16 Mx in the outer penumbra (r — 12 000 km), observed 
at a heliocentric angle of 8 = 30° and a spot angle of if/ = 0°. The flow speed inside the tube is Vo = lOkm/s which results in 
a line-of-sight velocity of vlos = vq ■ sin(0) = 5 km/s. Since the flux tube is aligned with the line-of-sight, the azimuth angle is 
zero. 



is looking at the center-side penumbra. For ip = 90° and 270°, 
the NCP vanishes since the line-of-sight is perpendicular to the 
flow. The azimuthal variation of the NCP is approximately an- 
tisymmetric with respect to the line connecting spot center and 
disk center. This behavior is characteristic of the infrared line 
Fe I 1564.8 nm and is known as the A(p effect (see Sect. l4.2l and 
M02 for details). 



4. Results 

4.1. NCP maps 

Figure [5] shows a comparison between NCP maps from obser- 
vations, the moving tube model, and a map that is constructed 
with the new VTUBE model. The upper row shows NCP maps 



of the Fe 1 1564.8 nm infrared line, the bottom row shows NCP 
maps of the Fe I 630.25 nm visible line: 

1. The first set of maps is constructed from measure- 
ments that have been performed at the German Vacuum 
Tower Telescope (VTT), feeding two polarimeters si- 
multaneously: the Teneri fe Infrared Polarimeter (TIP, 
iMartinez Pillet et all [l99<k for Fel 1564.8 n m and the 
Polar imetric Littrow Spectrometer (POLIS, beck et all 
120051) for Fe I 630.2 nm. The spot was at a heliocentric an- 
gl e of 8 - 30° . More details on the observations are given 
in lBecklll2005l) . 

2. The second set of NCP maps is based on the moving tube 
model, as described in detail in M02. 

3. The third set of maps demonstrates the capability of the new 
VTUBE model. For the atmospheric parameters, we use 
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Fig. 4. Stokes-V profi les and NCP for a horizon tal magnetic flux tube with a flux of T — 5 • 10 16 Mx, embedded in the penumbral 
atmospheric model of lBellot Rubio et al.N2006|) . observed at a heliocentric angle of 9 = 30°. The stack plot on the left side shows 
Stokes-V profiles for different positions along a radial cut at r = 12 000 km through the outer penumbra, specified by the spot 
angle if/. The right plot shows the corresponding NCP, i.e. the integral over the Stokes-V profiles. As expected, it is seen that the 
NCP vanishes for if/ = 90° and if/ = 270°, where the line-of-sight is perpendicular to the flow. 



results from two-component Stokes inversions. Analyses 
of such inversions have demonstrated that two interlaced 
components are ca pable of reproducing the spectropolari- 
metric data sets (c f. iBorrero et alj2006HBellot Rubio et alJ 
l2004lBeckl2005l) . The typical properties of these two com- 
ponents are summarized in Fig. [6] and used as input for 
VTUBE: From the inner to the outer penumbra, the mag- 
netic field strength of the tube component is about 500 G 
smaller than that of the background component and grad- 
ually approaches it as the radial distance from spot center 
increases. The flow speed increases from 6.5 to 8 km/s in 
the inner penumbra and stays constant in the outer penum- 
bra. The tube inclination gradually changes from being up- 
wards (a = 25°) to downwards (a = -10°), with a being 



measured relative to the horizontal. Assuming that these pa- 
rameters are axially symmetric, we employed the VTUBE 
model to calculate synthetic spectra for different viewing 
angles for the entire spot. For t he background atmosphere , 
we used the penumbra model o f lBellot Rubi o et alJ (|2006 ), 
and the tube was placed 50 km above optical depth unity 
with a magnetic flux of 5 x 10 16 Mx. The resulting NCP 
maps are shown in the right column of Fig. [5] 

The observations (left column) indicate that the NCP dis- 
tribution around the penumbra has two maxima and two min- 
ima in the infrared and is roughly antisymmetric with respect 
to the line-of-symmetry (line from spot center to disk center, 
indicated by a red arrow). The moving tube model (second col- 
umn from the left) is able to reproduce this fundamental dif- 
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Fig. 5. Maps of net circular polarization (NCP) for an heliocentric angle of 9 — 30°. Upper row: infrared line, bottom row: visible 
line. Columns from left to right: observations (VTT, Tenerife), NCP maps synthesized from the moving tube model, and maps 
from the new and generalized VTUBE model. The arrows point towards disk center. The absolute magnitude of the NCP is 
higher for the synthetic maps than for the observations since a filling factor of unity is assumed. The color coding for all maps 
represents the NCP, measured in pm, and the spatial scales are in Mm. 



ference. Since it is limited, however, to given snapshots of the 
MHD model, VTUBE serves to produce maps for other param- 
eter sets. In the third column we use a set of para meters as in - 
ferred from inversions of spectropolarimetric data jBeckl2005l) . 
The corresponding NCP maps reproduce the general symmetry 
properties, but have different radial dependences than the maps 



resulting from the moving tube model. When compared to ob- 
servations it is seen that the radial dependence on the limb-side 
is well reproduced, but that the center-side penumbrae differ in 
terms of radial dependence while the azimuthal dependence is 
comparable. In the case of the moving tube model, the sharp 
edges of the NCP contours in the inner penumbra are due to 
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Fig. 6. Radial dependence of flow and magnetic field of the two 
components. These values are used to construct the NCP maps 
in the right panel of Fig. [5] The left panel shows the radial vari- 
ation of the magnetic field strength of the strong background 
(dotted line) and weak tube (solid line) component. The mid- 
dle and right panels display the flow speed and the inclination 
of the tube component relative to the horizontal. 



flux tube embedded in a static background atmosphere. The 
NCP resulting from such an atmosphere with a flow channel 
can be caused by the gradient or jump in (a) the magnetic field 
strength (AB effect), (b) its inclination (Ay effect) and/or (c) 
the azimuth (A(p effect). As described in detail in paper M02, 
we showed that the A<p effect is the dominant one for the in- 
frared line Fe 1 1564.8 nm. In this case, the NCP is proportional 
to -\'los • sin(2A0) where vws is the line-of-sight velocity and 
A<p is the jump in the azimuth of the magnetic field between the 
flux tube and the background atmosphere. This term is plotted 
in the right panel and qualitatively reproduces the symmetry 
properties of the NCP displayed in the left panel. 

The second point we would like to make is that the ob- 
served azimuthal variations of the NCP in the infrared line al- 
ways show two maxima and two minima. Comparing this with 
models where the angle of the flux tube was varied, this finding 
rules out flux tubes that descend at an angle steeper than about 
a< -15°. 



the emerging flux tube which enters the photosphere and thus 
suddenly becomes visible in the two chosen spectral lines. For 
the VTUBE model calculations (right column), flux tube pa- 
rameters inferred from observations were used which indicate 
more gentle inclination angles of the flux tube component. This 
explains the smoother contours for the latter model. The abso- 
lute magnitude of the NCP is higher for the synthetic maps than 
for the observations since a filling factor of unity is assumed, 
i.e. each line of sight is passing throu gh the center of a mag- 
netic flux tube. A recent analysis by iBeckl d2005l) . based on 
Stokes inversions, indicates that the filling factor of the flow 
channel component is on the order of 40%. Since the NCP 
scales linearly with the filling factor, adopting this value sig- 
nificantly improves agreement between the magnitude of the 
observed NCP and the one c alculated wit h the new VTUBE 
model. However, the work of IBeckl d2005h also indicates that 
the filling factor varies throughout the penumbra, so a detailed 
analysis is needed in order to obtain a more qualitative compar- 
ison. In general, these results are very promising and have stim- 
ulated further modeling efforts. One open question is whether 
magnetic flux tubes dive back down into lower atmospheric 
layers in the outer penumbra and, if so, at which angle. 



4.2. The inclination of penumbral flux tubes 

Figure [7] shows the variation of the NCP along a circular cut 
in the outer penumbra for the infrared line Fe I 1564.8 nm. For 
tube inclination angles of a > -15°, the curve has two pro- 
nounced maxima and minima, while there is only one of each 
for smaller angles, i.e. more vertical downflows. 

There are two points to be made here. The first one 
is that the different shapes of the curves can be ex- 
plained by applying the simplified analytical model of 
lLandolfi & Landi DegFInnocentil ( ll996l) . In this model, a thin 
slab of moving plasma is located on top of a semi-infinite static 
atmosphere, and the two parts are permeated by different mag- 
netic fields. In M02 we showed that their findings can also be 
applied to the slightly more complex situation of a magnetic 



4.3. Effect of lower magnetic field strength in the flux 
tube 



Recent observation s by iBellot Rubio et all d2004l) and 
iBorrero et alJ d2004l) suggest that the magnetic field strength 
in penumbral flux tubes is lower than in its surroundings. One 
way of confirming this conjecture is to look for the signature 
of a discontinuity in the magnetic field strength in the NCP. 
Figure [S] shows the azimuthal variation of the NCP for both 
the infrared and the visib le line for the atmospheric model of 
IBellot Rubio et all d2006l) . The effect of lowering the magnetic 
field strength inside the flux tube is the following: For the 
infrared line Fe I 1564.8 nm, the NCP strongly increases in the 
center-side penumbra e [270°, 90°]). The dotted lines in 
the plots corresponds to a jump in magnetic field strength of 
AB = Bbg - B tl ,be = 200 G, the dashed lines corresponds to 
AB = 400 G. For the visible line Fel 630.25 nm, the effect is 
much smaller and consists of a small increase of the NCP in 
the center-side penumbra and a small decrease of the NCP in 
the limb-side penumbra. The explanation for this behavior is 
given in Fig. [5] The left panel shows the Stokes-V profiles of 
the infrared line, the right panel the Stokes-V profiles of the 
visible line (same line-styles as in Fig. [S}. Let us take a look 
at the center-side profiles around if/ = 0°: It can be seen that 
upon reduction of the field strength inside the tube, the two 
small humps of the flux tube component of the profile move 
closer together due to the reduced Zeeman effect. For the blue 
wing of the Stokes-V profile, this means that the small hump 
merges with the main lobe, which reduces the integrated area 
over the blue wing. For the red wing, on the contrary, the small 
hump becomes more and more disconnected from the main 
lobe. This increases the integrated area under the red wing, 
with the result that the total integral over the Stokes-V profile 
increases. In the case of the visible line, reducing the magnetic 
field strength increases the NCP just a little on the center-side 
(because the flux tube component of the Stokes-V profile is 
less split due to the Zeeman effect). On the limb-side, the NCP 
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Fig. 7. Variation of the NCP along a circular cut in the outer penumbra for the infrared line Fe I 1564.8 nm. For tube inclination 
angles of a > -15°, the curve has two pronounced maxima and minima, while there is only one of each for smaller angles, 
i.e. more vertical downflows. The different shapes of the curves can be approximated by the term -Vlos ■ sin(2A0) where vws 
is the line-of-sight velocity and A(f> is the jump in the azimuth of the magnetic field between the flux tube and the background 
atmosphere. 



is slightly reduced for the same reason (compare the profiles in 
the right panel of Fig. [5}. 

In conclusion, observing the limb-side NCP of sunspots in 
the Fe I 1564.8 nm line offers a promising way to identify a re- 
duced magnetic field strength in flow channels. If the NCP is 
significantly different from zero and does not show antisymme- 
try around the line connecting spot center and disk center, this 
is a strong hint for a lower magnetic field strength in the flux 
tubes. 



NCP (Fe I 1564.8 nm) 




180 
* [°] 



Fig. 8. Introducing a jump AB = B B q - B, u i, e in the magnetic 
field strength between the flux tube and the background in- 
creases the NCP. For the infrared line (upper panel), the effect 
is largest for the center-side penumbra. 



4.4. Effect of different background atmospheres 

In order to check whether the atmospheric background model 
has a strong influence on the NCP, we calculated the same 
model with different background atmospheres. The remaining 
parameters (r = 12 000km, a = 0°, 9 = 30°, AB = 0G) 
were kept constant. Figure[K)]shows stack plots of the resulting 
Stokes-V profiles and Fig.^Jdisplays the corresponding NCP 
curves. 

In the case of the infrared line Fe I 1564.8 nm (left panel), 
the profiles look similar in shape for all three atmospheric mod- 
els. 12iejin£_d££flis^e krgest for the JS94 model, followed by 
the lBellot Rubio et alJfcOOdl model and the lSchleicherl ( ll9761) 
model. This is simply a result of the temperature difference in 
the line-forming region between the three models: The JS94 
model is the coolest mod el and thus absorbs more strongly 
than the warmer mo del of [B ellot Rubio e t alJ ( 120061) and the 
even hotter model of ISchleicherT ( Il976l) . This is consequently 
reflected in the NCP curves which are similar in shape for all 
three atmospheres, but slightly larger in amplitude the cooler 
the model is. 

For the visible line Fe I 630.25 nm, the situation is slightly 
more involved. The right panel of Fig.lTolshows significant dif- 
ferences in the shapes of the profiles on the center-side. To be 
more precise, the unshifted components of the Stokes-V pro- 
files for the JS94 model differ from those of the other two mod- 
els in the sense that both wings are less pronounced than for 
the other two atmospheric models. The reason for this is that 
the JS94 model is a relatively simple theoretical one in which 
the temperature gradient almost vanishes above log(r) = -1, 
i.e. directly above the flux tube. For this reason, the line for- 
mation nearly ceases above the flux tube which results in the 
odd shapes of the unshifted Stokes-V components. A closer 
look at the line formation process reveals that the Stokes-V 
profiles for iff = 0° reach their final shape already around 
log(r) = -1.5 for this model, while this state is only reached 
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Fig. 9. The effect of a jump AB = Bsc - B lu i, e in the magnetic field strength between the flux tube and the background. The 
left panel shows the Stokes-V profiles of the infrared line, the right panel the Stokes-V profiles of the visible line. The solid line 
corresponds to AB = G, the dotted line to AB = 200 G and the dashed line to AB = 400 G. 



around log(r) = -2.5 a nd l og(r) - -3 for the models of 
iBellot Rubio et alJ feOOd and lSchleicherl(ll976l) . respectively. 

For the infrared line Fel 1564.8 nm, the vanishing temper- 
ature gradient of the JS94 model above the flux tube does not 
significantly affect the shape of the Stokes profiles since this 
line is formed in deeper layers of the atmosphere. 

Despite the fact that the Stokes-V profiles vary for the dif- 
ferent atmospheric models, the resulting NCP curves for the 
Fe I 630.25 nm line do not differ strongly. The slightly higher 
NCP for this model on the center-side penumbra can be at- 
tributed to the difference in line formation described above. 
From these findings, we conclude that the choice of the atmo- 
spheric background model can significantly influence the shape 
of the Stokes profiles, but does not change the global charac- 
teristics of the resulting NCP curves for the tested models. 
For comparison, Fig. H2lshows the observed azimuthal varia- 
tion of the NCP for the infrared line Fel 1564.8 nm and the 



visible line Fe I 630.25 nm. As noted above, the absolute mag- 
nitude of the observed NCP is lower than the one calculated 
from the synthetic profiles since the latter assume a filling fac- 
tor of unity, i.e. each line of sight is passing through the center 
of a magnetic flux tube. The exact calibration of the zero-level 
of the observed NCP is very difficult, therefore Fig. H2lshows 
the unshifted raw data. 



5. Summary 

We have presented a generalized geometrical model that em- 
beds an arbitrarily shaped flux tube in a stratified magnetized 
atmosphere. The new model is a versatile tool to calculate the 
spectral signature of flux tubes in the penumbra and especially 
to make predictions about the flow speed and tube inclination 
from observed maps of the net circular polarization. From the 
first applications of the new model, we find that 
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Fig. 10. The effect of different background atmospheres. The left panel shows the Stokes-V profiles of the infrared line, the 
right panel the Stokes- V profiles of the visible line. The solid li ne corresponds to the JS94 model, the dotted line to the 
iBellot Rubio et alJ J2006I) model and the dashed line to the model of lSchleicherl Jl976l) . 



(a) The inclination of downfiows in the outer penumbra must 
be shallower than approx. 15°. 

(b) Observing the limb-side NCP of sunspots in the 
Fe I 1564.8 nm line offers a promising way to identify a re- 
duced magnetic field strength in flow channels. 

(c) The choice of the background atmosphere can signifi- 
cantly influence the shape of the Stokes profiles, but does 
not change the global characteristics of the resulting NCP 
curves for the tested atmospheric models. 

In future work, we will compare synthetic NCP maps obtained 
from this forward-modeling approach with results from Stokes 
inversion maps of observational data in more detail. This will 
provide further insight into the structure of penumbral fila- 
ments and guide improvements in both modeling and inversion 
techniques. 
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Appendix A: Model description 

A.1. Defining the flux tube 

In order to construct a computationally efficient algorithm that 
is able to calculate the physical properties along the line-of- 
sight, we use the parametrized 3D model described below. In 
contrast to an explicit 3D model where one generates a three- 
dimensional box for all the physical properties and then inter- 
sects this box with a line-of-sight, our model generates the line- 
of-sight variables directly from geometrical considerations. 

We start with a two-dimensional definition of the flux tube 
axis, which is given by a polygon or spline function and is em- 
bedded in the (r, z)-plane as shown in Figs. IA.lt and b. This 
figure shows an example of a straight tube with a constant in- 
clination angle, a, but in principle, a is variable along the tube. 
A complete set of physical quantities relevant for the radiative 
transfer code is assigned to each point along the tube axis. The 
local radius of the tube, R tu be> is calculated from the condition 
that the magnetic flux, f, is constant along the tube, 



R 



tube 



(A.l) 



We then define the discretized LOS which intersects the tube 
at any desired point. The inclination 9 of the LOS is given by 
the heliocentric angle and its azimuth tfr corresponds to the spot 
angle. For all points L, (i = 1 , . . . , N tU be) along the LOS we 
calculate the nearest 3D distance, d^, to the tube axis as shown 
in Fig. IA.lt >. If d^ < R, u b e , the point L, lies inside the flux 
tube, otherwise it is associated with the penumbral background 
atmosphere. Based on this distance-dependent binary decision, 
we assign the physical properties to each point L, along the 
LOS. 

A.2. Assignment of the background properties 

Based on the coordinates (r,-,z,) of the points L, we can eas- 
ily interpolate the scalar quantities from tabulated background 
model values (temperature T, pressure p, density p, magnetic 
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Fig. A.l. Panel (a) shows a flux tube, inclined by an angle a and 
observed at a heliocentric angle 9 and a spot angle ip. Panel (b) 
shows the flux tube intersected by the LOS, projected into the 
(r, z)-plane. For each point L, on the LOS, the minimal distance 
di to the flux tube axis is calculated. If this distance is smaller 
than the local tube radius, the point L, lies inside the flux tube. 
Typically, the radiative transfer equation is integrated from an 
optical depth of log(r) = 2 out to log(r) = -5 (indicated as 
dashed lines for a vertical LOS). In panel (c), the relevant di- 
mensions for the interpolation of an arbitrary physical quantity 
X are depicted. First, a linear interpolation along the flux tube 
axis is carried out, thus leading to a value X ax j s (si) at distance 
Sj from Xj. Second, starting from the value X ax i S (sj) at the flux 
tube axis, one can prescribe an arbitrarily shaped radial depen- 
dence, which results in X Lj (sj, di). 



field strength B). In addition, one needs the magnetic field az- 
imuth <p an d its inclination y with respect to the LOS (see 
Fig. IA.2b and c). Let tws be the tangent vector along the LOS 



(directed into the spot) and B the magnetic field vector. The 
inclination y is then given by (compare Fig. IA.2t ) 

B ■ (-tws) , . ON 

r = arCC ° S lBU^- (A - 2) 
The minus sign arises due to the fact that we set y — 0° when B 
and tws are antiparallel, i.e. an inclination angle y e [0°,90°) 
indicates that the magnetic field points into the hemisphere of 
the observer. We assume positive magnetic polarity, i.e. the 
magnetic field vector is parallel to the flow vector inside the 
flux tube and points out of the solar surface for the magnetic 
background field. 



(a) 



(b) 




n I 


i 








(A.3) 



Fig. A.2. (a) Illustration of the local coordinate system (h, p), 
spanning the plane perpendicular to the line-of-sight. The az- 
imuthal angle <p with respect to the line-of-sight is shown in (b), 
whereas (c) shows the definition of the inclination between the 
line-of-sight and the magnetic field. 

In order to calculate the magnetic field azimuth <p, we have to 
project the magnetic field vector onto the plane perpendicular 
to the LOS as illustrated in Fig. IA.2h . We define p as the unit 
vector perpendicular to tws ar, d the z-axis (p is consequently 
located in the (x, y)-plane) and h as the unit vector normal on 
tws an d p- A closer look at Fig. IA.2b shows that <p is given by 

arctan (a/b) + tt/2 a > A b > 

- arctan (a/b) for a>0Ai<0 
arctan (a/b) + 2>n/2 a < Ab < ' 

- arctan (a/b) + n a < A b > 

where a = B ■ p and b - B ■ h. 

A.3. Assignment of the flux tube properties 

If a point L, lies inside the flux tube, the physical properties are 
interpolated from the values provided along the flux tube poly- 
gon. In particular, one first interpolates the quantity of interest 
(here symbolically termed X) linearly along the flux tube axis. 
If the two nearest flux tube neighbors of L, are the points j and 
7 + 1 with quantities Xj and Xj+\, respectively, we obtain for 
X ax i S in a distance s,- from j (compare Fig. IA.lb ) 

X axis (Si) = Xj + — ^- (X j+ 1 - Xj) . (A.4) 

Second, we have the freedom to analytically prescribe an arbi- 
trarily shaped radial dependence/ for X, 

X Li ( Si ,di) = f(X axis (si),di), (A.5) 
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with the limitation that f(X ax i S (si), 0) = X aX j S (Si). Such a de- 
pendence of the physical properties on the radius can e.g. be 
prescribed to model the cooling of a flux tube. In this study we 
do not investigate the impact of such radial dependences and 
use instead (as mentioned above) the thin flux tube approxima- 
tion, where the physical variables are radius-independent, 

f(X axis (s i ),dd=X axis (s i ). (A.6) 

The vector-valued quantities (i.e. the velocity and the magnetic 
field) are assumed to be oriented parallel to the flux tube axis 
and their strength is interpolated from the values given on the 
axis in the aforementioned manner. The projection onto the 
plane perpendicular to the LOS is carried out analogously to 
the previously described way (compare Fig. lA.2l and Eq. IA.3t . 

After all physical quantities have been projected onto the 
line-of-sight, they are passed to the radiative transfer code. 
Scanning the spot angle and the radial position of the line-of- 
sight within the penumbra, one can now easily obtain maps all 
four Stokes parameters or any derived quantities. 



